Rabi spectra - a simple tool for analyzing the limitations of RWA in modelling of the 
selective population transfer in many-level quantum systems 
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The controlled manipulation of quantum systems is becoming an important technological tool. Rabi 
oscillations are one of the simplest, yet quite useful mechanisms for achieving such manipulation. 
However, the validity of simple Rabi theory is limited and its upgrades are usually both technically 
involved, as well as limited to the simple two-level systems. The aim of this paper is to demonstrate a 
simple method for analyzing the applicability of the RWA Rabi solutions to the controlled population 
transfer in complex many- level quantum systems. 

PACS numbers: 3.65.Sq 



INTRODUCTION 



State selective manipulation of the discrete-level quan- 
tum systems such as atoms, molecules or quantum dots 
is the ultimate tool for many diverse fields such as laser 
control of chemical reactions, atom optics, high-precision 
metrology and quantum computing The elementary 
process on which many schemes of such manipulation 
rely are 'Rabi' population oscillations between the two 
selected levels of the system. In the lowest approxima- 
tion, these oscillations are induced by perturbing the sys- 
tem using coherent monochromatic external drive tuned 
into resonance with the desired transition. The name of 
this phenomenon stems from the name of the author of 
the first approximate analytic assessment of this process 



The original Rabi theory is simple to present and com- 
prehend, but has two drawbacks. The first is that due to 
the rotating wave approximation (RWA) , under which it 
is obtained, its validity is limited even in pure two-level 
systems. A number of papers has been published con- 
cerning the solution of population dynamics without the 
RWA which are all based on fairly involved math- 
ematical apparatus and hence quite un-intuitive. Also, as 
these methods are founded on perturbative approaches, 
the solution in the case of the extreme breakdown of RWA 
- which is extremely simple (although maybe not very 
useful) - is seldom mentioned. Secondly, the analytical 
approach to single-laser Rabi oscillations is generally con- 
strained to simple two level system, and hence its validity 
is necessarily limited when it is applied to more complex 
many-level systems. Since many (or, indeed, most...) in- 
teresting real systems actually are many-level systems, at 
least a rough assessment of influence of their structural 
complexity on the Rabi oscillations between the selected 
two levels should be done before applying them to com- 
plex systems. Finally, looking through the literature I 
have not managed to find a thorough discussion of the 



validity of RWA solution for the dynamics of a two level 
system driven by the perturbation whose frequency is far 
off resonance. 

In this paper all of these questions are addressed. In 
the first section, mathematically simple yet thorough 
analysis of the RWA limitations in the two-level systems 
is presented. In the second, a simple criterium (based 
on Rabi profiles) for the analysis of the validity of RWA- 
bascd solutions is established. This enables a quick es- 
timation of the validity of the Rabi solution for certain 
system parameters just by visually inspecting the shape 
of a simple lorenzian curve. Finally, the third section is 
devoted to discussion of how the concept of Rabi spec- 
trum may be used to estimate the validity of the two-level 
Rabi solution in a general many-level system. 

In all the calculations value of the Planck constant is 
taken to be h = 1 . 



II. ROTATING WAVE APPROXIMATION IN A 
TWO-LEVEL SYSTEM 

In this section a simple, yet thorough analysis of valid- 
ity of the Rabi RWA solution in the two level system is 
presented. The results obtained - majority of which are 
quite well known - will serve as strict mathematical basis 
for consequent discussions. 



A. Full dynamical equation 

Consider a two-level quantum system. Its dynamics is 
described by the time dependent Schrodinger equation: 



(flo(x) + y(i))V(x,t). 



(1) 



x denotes the set of all spatial coordinates of the system. 
Hq(x) is unperturbed time- independent Hamiltonian of 
the system, satisfying the eigenequation: 



'Electronic address: dbonacci@irb.hr 



H &(x) = E, &(x) ; (t = l,2), 



(2) 



2 



where (pi (x) and Ei are eigenfunctions and respective to- 
tal energies of the unperturbed system. V(x 7 t) is the 
homogeneous monochromatic perturbation: 



V(x,i) = CT(x)F cos(wi). 



(3) 



<r(x) is spatial operator odd with respect to central in- 
version (e.g. dipole, quadrupole, ...) through which the 
perturbation establishes coupling between the two levels 
of the system. Fq is the parameter determining pertur- 
bation intensity (which shall be referred to as the per- 
turbation intensity for short). u> is constant but tunable 
frequency of the external perturbation. Expressing the 
solution ip(t) as an expansion over eigenstates </f>i(x), 



i=l 



(4) 



Eq. transforms into a system of two coupled differ- 
ential equations for the time dependent expansion coeffi- 
cients Oj(t): 



-a 3 {t) = -zFocos^^a^I^e-'W' ; (j = 1,2) 

i—1 



Here Sij — 



(5) 

sji = sign(ujij), tOij = \u)i — u/j\ and Iij = 
J ^j(x)o'(x)^)j(x)dx. From the definition it is evident 
that I\2 — I21 and In — I22 = 0. Combining the two 
expansion coefficients into a two-component vector a(r) , 

expanding a cosine term as - — ^ and re-scaling the 

time variable to 



S12 — ^— t 



casts Eq. (JSJ) into convenient matrix form: 



d 



„2iTr 



,-2iAt 



-2iTr 



(G) 



a(r) (7) 



where A = ^^r 12 and T = 



Evolution of the population of a particular level is 
given by !!,(<) = Mr)| 2 . 



B. RWA solution and its limitations 

The rotating wave approximation hinges on the as- 
sumption that the complex exponential terms containing 
r in Eq. JJJ) have only minor (or even negligible) impact 
on the system dynamics and may hence be dropped in 
all further calculations. However, in order to fully ex- 
plore their influence, they shall be keep throughout this 
calculation. 

Now the following procedure is administered: first Eq. 
is differentiated with respect to r; then the same 

equation is used to eliminate da ~}^ that appears during 



the differentiation on the right-hand side; finally, the ma- 
trix equation is split into two ordinary coupled differen- 
tial equations for the coefficients ai(r) and 02 (t). The 
following expression is hence obtained: 



dr 2 



ai |2 (r) = 2a 1 , 2 (r)TA e T 2 %, 1 (r) 



(8) 



± Te 



±2iFt 



a 2 ,i(r) -2cos(2(r + A)r)ai, 2 (r) 



where the first index corresponds to the upper sign and 
the second to the lower. As no approximations have yet 
been introduced, this is still the full dynamical equation 
equivalent to Eq. Q. Initial conditions comprise the 
whole population of the system located in one of the two 
levels, with the other level being completely empty 

Keeping only the first term on the right-hand side of 
Eq. JHJ leads straight to the on-resonance Rabi RWA 
solution, sinusoidal population oscillations with period 
Oo = 7r and amplitude Aq = 1 (for details see e.g. |||). 
Whether or not the remaining terms on the right-hand 
side may be considered dynamically insignificant depends 
on the relation between the period ©0 and periods of os- 
cillating and complexly rotating parts of these remaining 
terms. If any of those is much smaller than Oo then the 
value of the elements on the right-hand side in which they 
are contained averages up to zero during any time inter- 
val small compared to the one in which ai(r) and 02(1") 
significantly change, and hence they play no dynamical 
role. Otherwise they have to be taken into consideration. 

Shifting the perturbation frequency just slightly away 
from the resonance, the second term in Eq. © enters 
the arena. Taken together, the first two terms lead to the 
'generalized' Rabi oscillations with period ©a = tt 



1 



and amplitude Aa = j^gz (for details see e.g. |(|)- So, a 
shift from the resonance increases the frequency of pop- 
ulation oscillations but at the same time reduces their 
amplitude. 

To proceed, note that T by its definition is always 
larger than A. Indeed, when perturbation is relatively 
close to the resonance - and if anywhere, the RWA so- 
lution should be valid in this spectral region - T >> A. 
The interplay of the relative dynamical significance of 
the remaining two terms on the right-hand side in Eq. 
© depends only on exact magnitude of T: if T >> 1, 
both of them rotate roughly at the same rate but as the 
ratio of their amplitudes is jj, the contribution of the 
first one overshadows the contribution of the second one; 
if however, the case is opposite, i.e. T << I , then the 
amplitude of the former is so much inferior to that of 
the later that the contribution of the third term can be 
neglected when compared to that of the fourth. 

Now consider the former case, that is T >> 1 so that 
additional dynamical features that generalized RWA so- 
lution fails to encompass are basically provided through 
the third term in Eq. JHJ. In full numerical solution 
for this case, these additional features are observed as 
the additional 'Bloch-Siegert' oscillations superimposed 
on the pure sinusoidal population transfer curve. Also, 



3 



an additional property of the full solution for this case 
which RWA Rabi solution cannot capture is the shift in 
the resonant frequency a; 12 • 

One way to quantify the limitations of the RWA Rabi 
solution for this case is to try to estimate what the am- 
plitude of these superimposed population oscillations is 
- as long as they are small enough, RWA solution can 
be considered as a decent approximation to the full solu- 
tion. This can be done as follows. The approximate non- 
RWA solution is sought as a sum aA (r) + ar (r) where 
aA (r) is the generalized Rabi RWA solution (taking care 
of the first two terms in Eq. (JSJ) and ar(r) the correc- 
tion due to the additional third term - remember that 
the fourth term is neglected because of the assumption 
that r >> 1. Inserting this sum into Eq. @ and as- 
suming that sla{t) >> ar(i") componentwise, the formal 
perturbative solution for the correction is obtained: 



a r (r) 



2iTr 



e 

-2iVr q 



aA(r)dr. 



(9) 



As r >> 1 holds, the frequency 8a of aA(r) will cer- 
tainly be much smaller than the frequency of complex 
exponentials inside the matrix. Hence, the obtained in- 
tegral can be evaluated between some To and r such that 
t — tq << 0a by taking almost constant sla(t) w aA (t"o ) 
outside the integral. The remaining integration is simple 
and finally yields: 



also inherently related to the inclusion of the 'non-RWA' 
contributions in matrix in Eq. Q. It should be kept 
in mind that in this case the perturbation is extremely 
strong. Under such circumstances the corresponding 
potential V(t) in dynamical equation (Eq. might 
need to be corrected with higher order interaction terms 
(quadrupole, octupole,...). Indeed, it may also happen 
that the structure of the original two-level system be- 
comes significantly distorted or even that the system 
breaks up completely. Hence this result might not be 
very useful from the point of controlled system manipu- 
lation. 

All in all, it may be concluded that the Rabi RWA so- 
lution is an excellent approximation to the exact solution 
of Eq.Q in complete spectral range < uj < 00 of per- 
turbation frequencies - and not just close to resonance 
UJ12 - as long as: 



r >> 1. 



III. RABI PROFILE ANALYSIS 



(12) 



In this section the notion of Rabi profile is introduced 
and it is demonstrated how the obtained condition for 
the validity of the Rabi RWA solution, T >> 1, can con- 
veniently be visualized in terms of it. 



a r (r) 







-2;rv 



aA(r ) (10) 



The frequency of the components of ar (r) is evidently 
much higher than that of the corresponding components 
of sla(t). As the amplitude of both matrix elements is 1, 
the ratio of amplitudes of the correction components and 
the generalized Rabi RWA solution components 



[ar i,2(t)| 
I q a 1,2(1-) I 



amounts p. Thus, while T >> 1, 8a(t) >> ar-(r) the 
RWA generalized Rabi solution is a good approximation 
to the full solution and the dynamical significance of the 
third term in Eq. JSJ may be neglected. 

The last remaining term in Eq. (JSJ has a fixed am- 
plitude equal to 4. It can be neglected if T + A >> 1, 
due to the zero-averaging effect of its rapid rotation, and 
since A can be reduced to zero, in the worst case this 
reduces to F » 1 - which is the same as for the third 
component. However, this term has a prominent role in 
the extreme situation, when T « 1. In such a case, 
the third term may be neglected due to its vanishing am- 
plitude, and since A < T even more so may the second 
term. Eq. JSJ then reduces to: 



-r^a 1)2 (r) 



2(l + coB(2(r + A)r))oi,a(r) (11 



which leads to the population oscillations with ampli- 
tude 1, but with temporally modulated period 0(t) ~ 
71 |i+cos(2(A+r)7-)| • ^ S w ith the third term, this one is 



A. Rabi profile of a single system transition 

Rabi profile of the transition between levels 1 and 2 
simply refers to the amplitude of population oscillations 
exhibited by the generalized Rabi RWA solution, con- 
sidered as a function of the driving frequency lo. It is 
parameterized with the total coupling strength (or cou- 
pling), D\2 = FqI\2 an d the resonant frequency of the 
system, OJ12: 



P12M 



1 



1 



{D12) 2 



(13) 



As can be seen from the Fig. 1, it is a simple bell- 
shaped curve centered on u)%2 with the half- width at half- 
maximum (or halfwidth) equal to Z?i2- 

Expressing the condition Ijl2|l in terms of the drive fre- 
quency, the resonant transition frequency and the cou- 
pling, a simple criterium for estimating the validity of 
the generalized Rabi RWA solution in the particular two- 
level system driven by the external homogenous coherent 
perturbation is obtained: D12 << 0J12 +lo. In the worst 
case this reduces to: 



D12 << W12 



(14) 



or, put in words, distance of the profile's center from the 
origin of the spectrum (u> = 0) must be much greater than 
the profile's half-width. We may also rephrase it as: 'At 
the origin of the spectrum the height of the profile must 
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FIG. 1: Rabi profiles of a single spectral line. In each graph, 
respective coupling D12 is indicated above the graph. Ob- 
serve how the profile's half-width increases with increasing 
coupling. The corresponding dynamics of a two-level system 
driven by the resonant perturbation is given in Fig. 2. On the 
abscissa, perturbation frequency uj is given in the units of the 
resonant frequency loyi of the system. 



be negligible'. If this condition is fulfilled, population of a 
two-level system driven with any desired frequency to will 
exhibit pure sinusoidal oscillations with amplitude equal 
to the height P\2{ui) of the Rabi profile at the selected 
frequency and period equal to: 

T12H = -^-Jp^J) (15) 

Fig. 2 exhibits four cases of the population dynamics 
in a resonantly driven two-level system. These are ob- 
tained as the 'exact' numerical solution to Eq. J2J (using 
the NDSolve method from MATHEMATICA software 
package). The four presented cases respectively corre- 
spond to Rabi profiles in Fig. 1, as indicated with the 
value of coupling D 12 above each graph. In the first case, 
the coupling is relatively weak so that the condition l|12|) 
is fully satisfied (actually, T = 50) and the pure sinu- 
soidal population oscillations are observed. The second 
case is a limiting one (r = 5): a minor rapidly oscillating 
Bloch-Siegert component may be observed superimposed 
on top of major sinusoidal oscillations. In the third case 
r = 1 so that the major discrepancy from the sinusoidal 
oscillations is expected, and is indeed observed, with the 
Bloch-Siegert and the Rabi components being roughly of 
the same order of magnitude. The last case (r — 0.1) cor- 
responds to lower limit of the extreme coupling (r << 1), 
and as predicted, the quasi-sinusoidal oscillations with 
time-dependent period are observed. 

Fig. 3 displays the numerical solutions for the popu- 
lation dynamics of a two-level system driven by the off- 
resonant perturbation. The Rabi profile (displayed in 
the top graph) is common for all four cases, with the 
coupling equal to D\2 = 0.05o;i2. By visual inspec- 
tion it may be checked that the condition for validity 
of the Rabi RWA solution is fulfilled, while analytically 
it amounts to T > 20. The four perturbation frequen- 
cies (w a = 1.050wi2, Wb = 0.850wi2, u c — 0.503wi2 
and ujd = 2.580wi2) are chosen so that Pi2(^a) = 0.5, 
P12K) = 0.1, Pia(w e ) = 0.01 and P u (u d ) = 0.001. 




FIG. 2: The exact numerical solution for the population dy- 
namics of the two-level system subjected to a resonant per- 
turbation. The quantity on the abscissa is time expressed 
through dimensionless variable tui2, and the ordinate is 
111,2 (t). Initially the complete population is located in level 
1. Populations Ili(t) and 112(f) are respectively represented 
with a full and a dotted line. Tick-marks are located at in- 
tegral multiples of the population transfer period, calculated 
from Eq. 1151 corresponding to the respective value of D12 . 



Vertical gridlines in bottom four graphs are drawn at in- 
tegral multiples of the predicted population oscillation 
time, determined from (|15|) . Horizontal gridlines are 
drawn at the predicted Rabi RWA value of the popula- 
tion oscillation amplitude determined from (|13[) . In the 
case a. the observed dynamics is clearly sinusoidal oscil- 
lations with the predicted amplitude and period. In the 
case b. the excellent agreement between the predicted 
Rabi RWA values and the numerical results is also clear. 
However, the minor Bloch-Siegert oscillations of ampli- 
tude y ~ 0.05 may be noted. In the case c. again the 
predictions are quite fair, although the amplitude of the 
Bloch-Siegert oscillations is now of the same order as the 
Rabi component so the dynamics is not purely sinusoidal, 
but is rather doubly periodic. In the last case, d. situ- 
ation is similar to the one in the case c, but as T is 
somewhat larger for this case, the agreement is some- 
what better. Altogether, the usefulness of the graphical 
criterium is clear. With narrower profiles situation is 
even better. 



B. Many- level system and the Rabi spectrum 

As a first instance of many-level system a three-level 
one shall be considered. It is actually the most important 
example, as all the consequent conclusions on the behav- 
ior of the many-level systems are a mere generalization 
of the findings from a three-level case. Let the internal 
structure of this system be such that the drive couples 
level 1 to level 2 with the total strength D12 = ^0^12 and 
level 2 to level 3 with the total strength D23 = F0I23, but 
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FIG. 3: Numerical solutions for the population dynamics of 
a two-level system subjected to non-resonant perturbations. 
On the top graph the common Rabi profile is displayed. Ver- 
tical lines indicate different perturbation frequencies for each 
of the four cases presented. Letters under the lines indicate 
which dynamics graph corresponds to each perturbation. On 
the four bottom graphs the corresponding population dynam- 
ics is displayed. Horizontal lines indicate the expected value 
of amplitude of population oscillations determined according 
to the corresponding height of the Rabi profile (Eq. 
in each particular case. Vertical lines mark the multiples of 
the oscillation period determined from Eq. 1)150. The excel- 
lent agreement of these two approximate expressions obtained 
within the RWA, and the 'actual' behavior represented by the 
numerical solution is quite clear. 



does not couple level 1 directly to level 3 (i.e. Irs = 0). 
The spectrum of the system hence contains two resonant 
frequencies, u)x% and W23. Initially let the whole popula- 
tion be located in level 2. The main goal is to transfer 
it to the level 1 minimizing both the transfer time and 
the loses to the level 3. Based on the results discussed in 
the previous section the dynamics of the three-level sys- 
tem might be predicted to some extent by inspecting the 
common plot of the two Rabi profiles P\2{u>) and P^siui) 
(top graphs in Fig. 4, 5 and 6). Such a plot containing 
all the Rabi profiles of the spectral lines which couple any 
of the two targeted levels (1 and 2 in this case) among 
themselves and to the rest of the system (just level 3 in 
this case) shall be referred to as the Rabi spectrum of the 
particular targeted transition. 

Assume that both profiles are 'narrow' so that for both 
of them Eq. (|14|l holds. Further assume that both pro- 
files are also 'narrower' then the distance between their 
centers, \uj\2 — ^23 >> Di2,D23- This situation is dis- 
played on top graph in Fig. 4, and it can always be 
achieved by sufficiently reducing the drive intensity, Fq. 
Now consider each of the two transitions separately. Con- 
sider tuning the perturbation to Wx2- As ^12(^12) = 1 
and ^23(^12) << 1, such perturbation would induce the 
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FIG. 4: Rabi spectrum and the numerical solution for the 
population dynamics of a three-level system. On the top 
graph the two Rabi profiles corresponding to the two reso- 
nant lines of the system are plotted alongside. The two spec- 
tral lines roughly correspond to the two of the low-laying ro- 
vibrational transitions in the HF molecule. The abscissa is ui 
(atomic units) and the ordinate is P12, 23(^0 (dimensionless) . 
Perturbation frequency is indicated with the vertical line on 
the Rabi spectrum plot, and it is set to the resonant frequency 
6l>12. The bottom graph presents the population dynamics of 
all three populations. The full line is Ili(t), the dashed line 
is Il2(t). The population of level 3 would be drawn as the 
dotted line, but it is so negligibly small at all times that it 
cannot be distinguished from the x-axis. 



complete population transfer oscillations on transition 
1 <-> 2 and almost negligible population oscillations on 
transition 2^3. Hence, if such a drive is applied to 
the three-level system, it might be expect that it would 
indeed produce complete population oscillations between 
levels 1 and 2 and would have almost negligible effect 
on level 3. Validity of this prediction can be clearly ob- 
served by inspecting the numerical solution presented in 
the bottom graph of Fig. 4. 

According to Eq. (|15|l . an increase in the speed of pop- 
ulation transfer can be achieved by increasing the cou- 
pling through increasing the perturbation intensity Fq. 
Suppose that system parameters (resonant frequencies 
and coupling integrals) are such that after the drive in- 
tensity is increased both profiles still satisfy Eq. H14fl . 
but that now the profile of the undesired transition 2-3 
broadens to such an extent that £23(^12) is still small but 
not negligible any more. This situation is displayed on 
the top graph of the Fig. 5. Considering again each tran- 
sition by itself, drive with frequency W12 would again pro- 
duce the complete periodic population transfer on tran- 
sition 1 <-> 2, but would also produce a clearly observable 
- albeit not complete - population oscillations on 2 <-> 3. 
So in the three level system we expect the population 
dynamics would change from that in the previous case as 
now both transitions 'react' significantly to the applied 
drive. Based on the discussion in the previous section, it 
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FIG. 5: Rabi spectrum and the numerical solution for the 
population dynamics of a three-level system. On the top 
graph the two Rabi profiles corresponding to the two reso- 
nant lines of the system are plotted alongside. The two spec- 
tral lines roughly correspond to the two of the low-laying ro- 
vibrational transitions in the HF molecule. The abscissa is ui 
(atomic units) and the ordinate is Pi2,23(w) (no units). The 
perturbation frequency is indicated with the vertical line on 
the Rabi spectrum plot, and it is set to the resonant frequency 
6l>12. The bottom graph presents the population dynamics of 
all three populations. The full line is Ili(i), the dashed line 
is Il2(t) and the dotted line is Tls(t). 
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FIG. 6: Rabi spectrum and the numerical solution for the 
population dynamics of a three-level system. On the top 
graph the two Rabi profiles corresponding to the two reso- 
nant lines of the system are plotted alongside. The two spec- 
tral lines roughly correspond to the two of the low-laying ro- 
vibrational transitions in the HF molecule. The abscissa is ui 
(atomic units) and the ordinate is Pi2,23(w) (no units). The 
perturbation frequency is indicated with the vertical line on 
the Rabi spectrum plot, and it is set to the resonant frequency 
W12. The bottom graph presents the population dynamics of 
all three populations. The full line is Ili(i), the dashed line 
is Fl2(t) and the dotted line is II3 (t) . 



can be suggested that it would roughly comprise of most 
population being exchanged along transition 1 <-> 2 with 
the remaining fraction being now also exchanged among 
levels 2 and 3. This is clearly observed from the numer- 
ical solution displayed on the bottom graph of Fig. 5. 
It can also be observed that the maximum amplitude of 
the perturbing level's population equals almost exactly 

By further increasing the perturbation intensity (Fig. 
6), a simple sinusoidal pattern of oscillations on the tran- 
sition 1 <->• 2 disappears completely. Hence the perturba- 
tion this strong cannot be employed to drive the efficient 
population transfer along the transition 1 «-> 2. This 
suggests that in a general many-level quantum system 
there is an upper limit to speed of the efficient popula- 
tion transfer (i.e. such that the great majority of the 
population is transferred) along each particular transi- 
tion. This limit is determined exclusively by the internal 
properties of the system. Consider two levels a and (3 of 
a many-level system that are directly coupled by the ap- 
plied drive, i.e. such that I a p ^ 0. The efficient transfer 
along the transition a <-> f3 is feasible as long as none of 
the profiles from that particular transition's Rabi spec- 
trum has significant height at the resonant frequency of 
that transition - except, of course, the Rabi profile of 
that selected transition (which is 1 at its own resonant 
frequency). Certainly, it should also be kept in mind 
that the Rabi profile of that selected transition has to 
have negligible height at the origin, i.e. it has to satisfy 



the condition (|12|l for the validity of pure RWA solution 
at that transition. 



IV. CONCLUSION 

The most important result of this paper is the estab- 
lishment of the simple graphical criterium for the validity 
of the approximate two-level RWA solution in modelling 
the population transfer in a general many-level discrete 
quantum mechanical system. The first hint of Rabi spec- 
trum method that I am aware of appeared quite some 
time ago - Fig. 5 in that paper is actually the Rabi 
spectrum plot. Looking at Fig. 6b in the same paper, 
minor oscillations can be noticed superimposed on top of 
almost pure Rabi oscillations of one of the populations. 
These minor oscillations gave me the idea that Rabi pro- 
file plots might be used as an indicator as to when the 
pure Rabi oscillations, which are a straightforward two- 
level RWA solution, fail to provide the satisfactory solu- 
tion to the complete dynamics of a more complex many- 
level system. 

Although the presented analysis is mathematically 
quite simple, it nevertheless enables a very good order- 
of-magnitude estimates for maximum driving field inten- 
sities that can be used to efficiently drive the selected 
transition. Rabi spectra method hence hints on one in- 
teresting application. Namely, knowing the details of the 
internal structure of a particular quantum system, the 
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Rabi spectrum for each transition in that system can 
be calculated. This data can then be used to quickly 
quantitatively estimate the lower limit for clean popula- 
tion transfer time that can be achieved on any particular 
transition. Finally, a numerical scheme can be build to 
search for the quickest pathway through the system along 
which population can be transferred between any two se- 
lected levels, and lower limit of the time needed can be 
estimated. 

Most of all, Rabi spectra method offers conceptual sim- 
plicity and intuitiveness which are both of utmost impor- 
tance when whole subject is presented either to interested 
non-experts or to newcomers in the field. Hence it can 
be considered very suitable both as a quick lab-reference 
tool as well as a simple didactic introduction into much 



more involved strict quantitative theory. 

The full analytical treatment of the selective popula- 
tion transfer mechanism provides many additional inter- 
esting results. It can thus be shown Q that by using 
analytically determined chirp (i.e. time dependent vari- 
ation of the perturbation frequency) instead of resonant 
frequency to drive any particular transition, more effi- 
cient population transfer can be achieved. 
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